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ABSTRACT 

We examine in detail several aspects of the physics of accretion disks of 
possible relevance to the outburst mechanism of the black hole X-ray transients. 
We adopt the one dimensional, time dependent model described in detail 
by Cannizzo, Chen, & Livio with parameters appropriate for a system such 
as A0620-00. We investigate (1) the effect of the grid spacing, utilizing a 
logarithmic radial spacing Ar oc r in addition to the spacing Ar oc ^/r^ (2) 
the dependence of the local flow speed of gas within the hot part of the disk 
on radius and time during the time of the cooling wave propagation, (3) the 
shape of the outburst light curve as a function of the triggering location for 
the instability, (4) the long term light curves of outbursts taken from trials in 
which complete cycles of quiescence and outburst are followed, both with and 
without including the effect of evaporation or removal of matter from the inner 
edge of the disk, and (5) the strength of the self-irradiation of the outer parts 
of the disk by the X-rays from the inner disk. Our primary findings in each of 
these areas are that (1) low resolution runs taking A?" ~ 20 grid points using 
the logarithmic spacing produce decay time scales which are artificially slow by 
factors of ~ 2 — 3 and slower-than-exponential, (2) the deviation from steady 
state within the outer part of the inner hot disk appears to be in accord with 
the discussion given in Vishniac & Wheeler — far from the transition front the 
fiow speed is ~ aCs{h/r), whereas at the interface between the transition front 
and the cold disk the fiow speed is ~ ac^, (3) the outburst triggering location 
must be > lOOrjnner for the rise time of the resulting outburst to be as short 
as is observed in the standard, bright systems, (4) the long term light curves 
using the standard model produce frequent outbursts which are triggered near 
the inner disk edge and which have slow rise times, and the long term light 
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curves calculated assuming evaporation of matter from the inner disk exhibit 
outbursts with longer recurrence times and somewhat (but not significantly) 
shorter rise times, and (5) for a system with parameters relevant to A0620-00, 
our "standard" system, irradiation is not a dynamically significant effect — in 
accord with recent results by van Paradijs. 

Subject headings: accretion, accretion disks — binaries: close — black hole physics — 
stars: individual (A0620-00) — X-rays: bursts — X-rays: stars 
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1. INTRODUCTION 



The black hole X-ray binaries (BHXBs) are interacting binary systems in which a 
Roche lobe filling K— or M— type secondary transfers matter through the inner Lagrangian 
point into orbit around a black hole accretor. The systems are therefore similar in many 
respects to the cataclysmic variables, the difference being that the accretor is not a white 
dwarf. (For a critical comparison of soft X-ray transients and dwarf novae, see van Paradijs 
& Verbunt 1984.) The BHXBs were discovered because of large outbursts seen first in 
X-rays which lasted weeks to months. The systems were subsequently studied optically, 
and found to be binary systems with orbital periods of several hours to several days, and 
with large mass functions /(M) > 3Mq implying black hole primaries (van Paradijs & 
McClintock 1995). 

The accretion disk limit cycle mechanism which was originally put forth to account for 
the dwarf nova outbursts (Meyer & Meyer- Hofmeister 1981; for recent reviews see Cannizzo 
1993a, 1998b; Osaki 1996) should work in other interacting binary stars. The model 
basically operates due to an inherent hysteretic relation between the vertically integrated 
viscous stress i^E and the surface density E at a fixed radius in the disk. The quantity u is 
the kinematic viscosity coefficient u = 2aP/{3flp), and a is the Shakura & Sunyaev (1973) 
parameter characterizing the magnitude of the viscous stress. When equilibrium solutions 
are plotted as loguT, (or logTeflfective) versus logE, one sees an 5"— curve relation. The lower 
and upper branches of the S which have positive slope are thermally and viscously stable, 
whereas the middle negatively sloping portion of the S is unstable. During quiescence, the 
gas in the accretion disk is far from steady state: matter accumulates at large radii in the 
disk and almost none reaches the central object. The gas is neutral, and the temperature 
and surface density values it possesses place it along the lower stable portion of the S curve. 
As matter continues to pile up, at some point a critical surface density Emax corresponding 
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to the local maximum in the S curve is exceeded, and a heating instability is initiated 
which ultimately takes all or most of the disk to the hot state wherein the gas is ionized 
and resides along the upper stable branch of the S curve. The rate of viscous evolution is 
much greater than before, and the E(r) profile shifts to a near steady state configuration 
in which E(r) decreases with radius, in contrast to the quiescent configuration in which 
S(r) increases with radius. After some material has accreted onto the central object and 
produced an outburst, the surface density decreases until at some point the local value drops 
below another critical value T^uiin characterizing the lower bend in the S curve. This always 
happens near the outer disk edge, and initiates a cooling wave which propagates inward to 
smaller radii, shutting off the flow and reverting the disk material back into neutral gas. If 
this model is correct, one cannot avoid having the limit cycle instability operating — unless 
the mass transfer rate feeding into the outer edge exceeds a critical value (Smak 1983), 
or irradiation is able to keep the entire disk permanently ionized (Tuchman, Mineshige, & 
Wheeler 1990; van Paradijs 1996). Not long after its introduction, various workers began 
to apply the model to the soft X-ray transients — a class which includes systems with both 
neutron star and black hole accretors (Cannizzo, Wheeler, & Ghosh 1982, 1985; Lin & 
Taam 1984; Huang & Wheeler 1989, hereafter HW; Mineshige & Wheeler 1989, hereafter 
MW). 

Cannizzo, Chen, & Livio (1995, hereafter CCL) investigated in detail the decay of the 
hght curves accompanying the outbursts in the BHXBs. They discovered that, if a cooling 
wave is responsible for the decay of the outbursts, then the observed exponential decay 
of the soft X-ray fluxes in the BHXBs constrains the viscosity parameter a to have the 
form a — Q;o(/i/r)^'^, where ckq — 50 if the primary (i.e., black hole) mass in the BHXBs is 
~ 1OM0. This form was introduced by Meyer & Meyer- Hofmeister (1983, 1984). The value 
of ~ 50 for the proportionality constant is mandated by the observed ~ 30 d e— folding 
decay times for outbursts in the BHXBs. The constraint is actually on Mi/aQ. A recent 
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compilation of inferred black hole masses for BHXBs by Bailyn et al. (1996) seems to 
indicate that the distribution is strongly peaked at ~ 7Mq, whereas the value ao ~ 50 
was based on assuming Mi ~ lOM©. If the work of Bailyn et al. is borne out, then we 
would need ckq ~ 35. To maintain continuity with CCL, however, we keep Mi = lOM© and 
CKo = 50. 

Vishniac & Wheeler (1996, hereafter VW; see also Vishniac 1997, hereafter V97) 
present an important study which provides a framework for understanding CCL's findings 
based on considerations of the deviation from steady state flow conditions within the hot, 
inner part of the accretion disk. Their key insight is that the flow within the hot part of 
the accretion disk which exists at small radii during the decay from maximum light after 
an instability has been triggered exhibits a deviation from steady state conditions as one 
approaches the interface with the cold portion of the disk. Their approach to the problem 
represents a fundamental overthrow of the old idea of expressing the transition front speeds 
in terms of purely local conditions (Meyer 1984). Part of our motivation in this work is to 
understand and to test their model. 

CCL's model for the limit cycle instability had an advantage over previous studies 
in that their numerical resolution was much flner than in other models, and they could, 
through systematic testing, place limits on the extent to which certain aspects of the model 
could be trusted. In this work we continue with this line of inquiry by running a series of 
models based on the numerical code described in previous works (Cannizzo 1993b, 1994; 
hereafter C93b and C94, respectively; CCL) to investigate various aspects of the model 
which were studied by CCL. In particular, we study (i) the effects of grid spacing, (ii) 
the analytical model presented by VW, (iii) the shape of the outburst versus instabihty 
triggering radius, (iv) long term light curves covering quiescence and outburst, both with 
and without incorporating evaporation of gas from the inner edge of the disk, (v) and the 
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strength of self-irradiation of the outer disk by flux from the inner disk. In the Discussion 
section we attempt to synthesize the lessons that were learned in the previous sections 
into (somewhat of) a coherent picture of the operation of the limit cycle instability in the 
BHXBs. 

2. MODEL CALCULATIONS 

For our standard model we use that presented in CCL. The code has been described 
extensively (e.g., C93b, C94, CCL): it solves explicitly for the time dependent evolution of 
surface density S(r, i) and midplane temperature Tj^idii^jt) by simultaneously calculating 
small corrections in E and T^id based on the equations of mass and energy conservation. 
The scalings parametrizing the steady state equilibrium disk structure arc given in C93b 
and CCL, and are based on Cannizzo & Wheeler (1984), Meyer & Meyer-Hofmeister (1981), 
and Pringle, Verbunt, & Wade (1986). 

C93b carried out a series of tests using the time dependent model to determine which 
terms are dominant and which are less important in their effect on the long term light 
curves. They also looked at the effect of varying the number of grid points. The results 
of the extensive testing done by C93b have been generally confirmed in a recent work 
which uses a more advanced time dependent disk instability code that solves not only the 
vertically averaged viscous and thermal energy equations but also the radial Navier-Stokes 
equation and therefore includes radial pressure gradients and departures from Keplerian 
flow within the transition fronts (Ludwig & Meyer 1997). 

We assume a central object mass Mi = IOMq, a feeding rate into the outer disk 
3 X 10~^^Mq yr~^ (McChntock et al. 1983), a viscosity parameter a — 50(/i/r)^-^, an inner 
radius rinner = 10^ cm, and an outer radius router = 10^^ cm. We use a grid spacing Ar oc y/r 
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in sections 2.2, 2.3, and 2.4, and a spacing Ar oc r in sections 2.1 and 2.5. In the three 
sections for which Ar oc -Jr, we take = 1000 grid points where we look at the properties 
of individual outbursts, whereas we take N = 500 grid points for the sections which treat 
complete cycles of outburst and quiescence. For section 2.1, we take N — 2"*(10) + 1, where 
1 < m < 5, and for section 2.5 we take m — 3 (i.e., N — 81). 



2.1. Logarithmic Grid Spacing 

CCL utilized a grid spacing Ar oc ^/r. They discuss earlier time dependent work by 
MW, Ichikawa, Mineshige, & Kato (1994; hereafter IMK), and Kim et al. (1994, 1996; 
hereafter K94 and K96, respectively) in which a logarithmic spacing Ar oc r is adopted, 
taking ~ 20 — 40 grid points. Some of the discussion in CCL concerning the number of 
grid points required to get believable results was speculative in the sense that CCL did 
not attempt to run models with the logarithmic spacing to substantiate their criticisms, 
but rather they extrapolated from what they had found using their high resolution trials. 
It seems worthwhile to look briefly at this issue of the logarithmic grid spacing so as to 
quantify the requirements for good results. Figure 1 of K94 shows a trial light curve 
using router /'"inner = 10'', Mi = 3.2Mq, and a = ao{h/ry-^. To compare our results more 
directly with K94, we have generalized our code so that the logarithmic spacing can be used 
(Cannizzo 1996b). 

Figure 1 shows the results of a series of decays for N = 21, 41, 81, 161, 321, and 501 
grid points using the logarithmic spacing. For ease of viewing the last run was started from 
different initial conditions so as to be offset slightly from the first five. Several aspects of 
the decay are apparent. The lower resolution runs have a decay which is too slow. For the 
N — 21 trial, the decay constant ~ 60 — 100 d, whereas in the limit of large we see 
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te — >~ 30 d. Note also that in the low resolution runs one sees ghtches in the light curves 
and te curves associated with the turn-off of individual grid points ( "turn-off" meaning the 
initiation of the cooling transition). Furthermore, te increases with time within a given run 
for N small, whereas we know from CCL and VW that for the parameters used, namely 
a oc {h/ry-^, the decay constant tg should not vary significantly in time. Therefore, the 
detailed decay rate in studies such as those mentioned earlier which used N ^ 20 grid points 
with the logarithmic spacing have decays which are spuriously slow by a factor ~ 2 — 3, 
and for which the functional form for the decay is more concave upwards (i.e., tg increasing 
with time) than it should otherwise be. 



CCL present numerous tests of the limit cycle model using a — ao{h/r)'^ for the 
Shakura-Sunyaev a parameter. They show the results of experiments for different values 
of tto and n and find that for runs with n close to 1.5, the decays of the outbursts have 
a nearly exponential form, as required by observations. CCL attempt to understand this 
result by casting the speed of the coohng front in terms of its width. VW show that this 
viewpoint is not valid: the width of the cooling front is an artifact of its speed, and not the 
cause. A key part of VW's insight into the physics of the cooling front and outflow from 
the hot, viscously evolving part of the disk is that, far away from the front the local flow 
speed Vr ~ acs{h/r), whereas in the front the flow speed Vr ~ aCg. VW show how this leads 
to a front speed 



where the subscript / denotes values within the cooling front, 1/q = 

1 -|- (1 + n/2){ai/[l — nbi/2]), n is the exponent associated with the a scaling, and 



2.2. A Close Look at the Cooling Front 




(1) 
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Oi and 61 are logarithmic derivatives characterizing the steady state scahng associated with 
the ionized state of the disk, ai = 0? log T/rf log S and bi = dlogT/dloga (VW; CCL). 
In other words, the scaling for the midplane temperature associated with the hot disk 
is Tjaid E"iQ;''ir~'^^ For the usual Kramer's law opacity, ai — ci — 3/7 and 61 = 1/7 
(C93b; CCL; VW). The midplane temperature associated with the local minimum in the 
logTgfj — logS steady state curves Tmm oc a~^'^^'^il~^^'^^ (VW). Assuming a oc (h/r)"' and 
using the fact that the sound speed in the front c/ oc ^/T^ yields dlogT^m/ dlogr — 
[9/140 - (l.l/7)(n/2)]/[l + (l.l/7)(n/2)]. (This last expression which is eqn. [34] of VW 
corrects an error in eqn. [4] of CCL, namely the r exponent should have positive rather 
than negative sign. In addition, in eqn. [4] of CCL there should not be a "K" on the right 
hand side.) 

There are several predictions of VWs theory which are amenable to testing with our 
time dependent code. We examine first the evolution of Vr/ aCg in space and time. Figure 2 
shows the variation of the local values of Vr/otCs and h/r with radius for eight time steps 
spaced 10 d apart during the early decay of a model outburst. The local flow speed is in 
fact ~ aCs{h/r), the usual viscous flow speed, some distance upstream from the transition 
front. VW speculate that this point occurs only slightly interior to r/, whereas we find it 
happens at ~ (l/2)rj. VWs derivation is predicated on the fractional disk radius at which 
this point occurs remaining constant during the decay, which it does. Also, the flow speed 
does not attain ~ acg at the inner edge of the front, as posited by VW, but rather at the 
outer edge. At the inner front edge the flow has only attained a speed ~ {l/6)acs. (In 
looking carefully at Fig. 2 one might have imagined a scaling with \Jh/r in the value of 
Vr/cs at the inner edge of the front, but an extension of the run used to generate Fig. 2 
in which the cooling front is followed down to ~ 10^ cm shows that the shoulder in the Vr 
curves which lies at the inner edge of the cooling front maintains a constant level in terms 
of Vr/cg.) VW note that CCL's results hinted that ~ {l/6)aCs at the inner front edge, and 
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the more detailed model presented by V97 uses this information as a refinement. Figure 3 
shows, for the eight time steps indicated in Fig. 2, the values of the radii at which the local 
flow speed goes to zero and the local flow speed equals the standard viscous value acs{h/r), 
both given in units of the radius associated with the outermost grid point in the hot state 
'"(^hot)- In the early stages of the outburst the ratio associated with the Vr — radius is 
> 0.6, but later it approaches a limiting value ~ 0.4. The formalism presented in VW 
which assumes Vr ~ acs for the flow speed at the inner edge of the cooling front flnds the 
Vr — fractional radius to be ~ 0.23, whereas the formahsm in V97 which is "corrected" 
by assuming Vr ~ {l/6)acs at the inner edge of the cooling front finds the Vr = fractional 
radius to be ~ 0.36, in accord with our results. The ratio associated with the point at 
which viscous fiow is reached asymptotes to ~ 0.55, smaller than that mentioned in VW. 
The deviations seen in the first few time steps from the limiting values are an artifact of 
the transitory phase near outburst maximum when the outwardly moving heating front 
has stopped near the outer disk edge and begun to reflect inward as a cooling wave. This 
artifact is also seen in the width of the transition front, which is quite broad near the end 
of the peak of the outburst just as the decay is beginning (see the upper panel in Fig. 5 of 
Cannizzo 1996c). Last, V97 note that the rate of mass outflow at the front is about three 
times the rate of accretion onto the central object. This is in agreement with Fig. 9 of 
CCL, and also Fig. 2c of C94. In summary, the discrepancies we find with VW are minor 
and it seems they have captured the essential physics. 

A further prediction of VW is contained within their expression for the cooling front 
speed. Their formula is more complicated than one might have guessed from the discussion 
in CCL. One sees for instance a dependence on oi. From equation (1) we have that 
Vf oc c'j'^'^'^^r^^'^"'^/'^ _ For the standard model parameters given in CCL and reiterated 
earlier, the optimal n value for which Vf oc r/ is not 1.5 but rather 1.656. In Figure 4 we 
show the results of varying n in much smaller increments than was attempted in CCL. 
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Eleven curves are given corresponding to n — 1.50, 1.52, 1.70. Rather than guessing 
an initial configuration close to outburst maximum and then following the decay, as was 
done in CCL, here we initiate an outburst at a small radius and wait until the disk has 
gone nearly entirely to the high state and then begun its decay. In this manner the disk 
adjusts to its desired state as the outburst decay starts, and we are not left with a long 
transient near the beginning of the decay as in CCL. In Figure 4 we see that, contrary to 
the discussion given in CCL, n — 1.5 does not yield a precisely exponential decay. Our 
optimal value is about 1.63, which is close to the value 1.656 given by VW's analytical 
formalism. Note the small dynamic range in e— folding times represented along the axis. 
The observational uncertainties and intrinsic source fluctuations would probably not allow 
one to distinguish among the range of n values shown here, so this result is not of direct 
relevance observationally. The motivation is solely to confirm VW. It would be a worthwhile 
exercise to analyze the decays of dwarf novae and X-ray novae in a consistent fashion to see 
how precise of a constraint one may be able to place on n, or, more broadly, how essential 
is the functional form a — ao{h/r)'". 

For completeness we show in Figure 5 the variation of the width of the front during 
the rise and decay of an outburst. The outburst begins at a small radius, so we see first 
the outward propagation of the heating front, then the inward movement of the cooling 
front. As noted by CCL, the width of the cooling front can be approximated by w ~ Vhr 
or w/r ~ \Jh/r. VW show that both the coohng front width and speed of propagation 
are set by the departure from steady state discussed earlier and are not causally linked. 
It is interesting that the fractional width of the outward moving heating front is given 
phenomenologically by w/r ~ {h/r)^^'^^. Ludwig & Meyer (1997) find somewhat broader 
heating fronts. Their Figure 14 is interesting because it plots the heating and cooling front 
velocities versus radius on a common scale. VW put forth an analytical theory for the 
period of decaying light associated with the inward moving cooling front; no such theory 
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currently exists for the period of rising light associated with the outward moving heating 
front. The sharp spike in E accompanying the heating waves is caused by the high viscosity 
material pushing into the region of lower viscosity material (Lin, Papaloizou, & Faulkner 
1985). The fact that it is narrower than the coohng front must be an artifact of this 
pushing, but further work will be needed to determine a valid numerical parametrization of 
the speed in order to lead ultimately to a successful theory of the heating front. 

2.3. Outburst Light Curve Versus Triggering Location 

It has been known for some time that, within the context of the limit cycle model, 
outbursts triggered near the inner edge of the disk tend to show slow rise times, while 
outbursts triggered at somewhat larger radii generate faster rise times (Smak 1984; 
Cannizzo, Wheeler, & Polidan 1986, hereafter CWP). The reason for the difference in rise 
times is detailed in CWP. For most of the bright BHXB outbursts, one sees fast rise times 
(of order a few days) in outbursts lasting ~ 100 — 200 d. It is therefore of interest to find out 
what is the minimum critical radius such that a fast rise will be obtained. Figure 6 shows a 
series of outbursts generated by starting with a smooth S(r) distribution resembling that 
expected at the end of the quiescent (i.e., accumulation) phase, and adding to this profile 
a small spike whose maximum exceeds the local value of Emax- This is carried out by 
superposing two functions, as shown in the top panel of Fig. 6. The functions are 



where Tmid =10^ cm is the geometric average of the inner and outer radii, reenter is the 
central radius of the narrow Gaussian used to trigger the outburst, w — rcenter/3, and 





and 




(3) 
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is the local maximum in E from vertical structure computations (Cannizzo & Wheeler 
1984). For a given trial we prescribe the initial surface density distribution by taking 
E(r) = max(5'broad, 'S'narrow)- Fivc trials were run, corresponding to log reenter ( cm) = 8, 9, 10, 
10.25, and 10.5. (For ease of viewing we do not show the 10.25 and 10.5 curves in the top 
panel of Fig. 6.) Note the stark difference in the rise time properties between the first two 
and last three curves. The rise is slow for the log reenter (cm) = 8 and 9 curves, whereas it is 
fast for the log reenter (cm) = 10, 10.25, and 10.5 curves (Cannizzo 1996b). The instability is 
triggered at t — for all five trials, but due to the finite travel time of the heating front, the 
effect of increased viscosity and mass flow rate is not felt at the inner disk until ~ 30 — 50 d 
later for the three trials with triggering at large radii, whereas the effect is felt much sooner 
for the trials triggered at small radii. The difference in the heating wave properties is shown 
in Figure 7. The two panels show the evolution of surface density, in 5 d steps, for the 
logr(cm)= 9 and 10 runs. These two span the transition region of parameter space between 
slow and fast rise. For the outbursts triggered at large radii, the amplitude of the outburst 
is larger because more mass is contained within the wave of surface density enhancement 
associated with the heating front. The slow rise times associated with the logr(cm)= 8 
and 9 curves are caused by the fact that, to enhance appreciably the rate of mass flow 
at the inner disk edge, the outwardly traveling heating front must access the bulk of the 
matter which is stored in the disk at large radii, between the Lubow & Shu (1975) radius 
(corresponding to that point in the disk where the specific angular momentum associated 
with the gas leaving the LI point at the surface of the secondary star equals the Keplerian 
value in the disk) and the outer disk edge, and this matter, after having made the transition 
from neutral to ionized gas, then has to viscously evolve (both to smaller and larger radii). 
The physics associated with this flow is discussed in CWP and C93a. 
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2.4. Long Term Light Curves 



As noted by many authors, the standard Umit cycle model fails to account for 
observations of systems in quiescence (see Lasota 1996a, 1996b). One observes a significant 
EUV and X-ray flux which we infer to be the result of accretion during quiescence, and 
yet if we are to take the theory literally it predicts a quiescent accretion rate which is 
negligibly small. The small accretion rate is mandated by the requirement that E < Emax 
and Teff < ^eff (^max) for all radii, in particular r ~ Tinner- McClintock, Horne, & Remillard 
(1995) found evidence for low level accretion ~ 10^° — 10^^ g (if the efficiency of accretion 
e ~ 0.1) in the quiescent state of A0620-00, but this is still much greater than would be 
allowed in the "pure" version of the standard accretion disk model in which the inner disk 
edge extends all the way to the last stable orbit (Lasota, Narayan, & Yi 1996). The critical 
accretion rate associated with E^in is M(Eniin) ^ 10^ g s~^ (nnner/lO^ cm)^'^ for BHXBs 
(CCL). There is nothing in the model, however, which requires Tinner = 6GMi/c^. Time 
dependent computations of the accretion disk limit cycle instability have shown that the 
model still works well if the inner edge is truncated (e.g., Angelini & Vcrbunt 1989; C93b). 
There may exist some physical mechanism extrinsic to the standard limit cycle model 
which serves to evaporate or take material from the inner disk and transfer it more or less 
directly onto the accretor (e.g., Meyer 1990; Meyer & Meyer-Hofmeister 1994; Liu, Meyer, 
& Meyer-Hofmeister 1995). Such a mechanism may or may not fit in with the general 
paradigm of accretion onto black holes (Chakrabarti 1996a, 1996b; Narayan, McClintock, 
& Yi 1996). Workers are beginning to consider the effect such a process would have on the 
limit cycle model in terms of how the altered surface density distribution in quiescence — in 
particular, the lack of an inner disk — would change the properties of the outburst (Lasota 
et al. 1996; Hameury et al. 1997, hereafter HLMN). 

When we run complete cycles using the a — 50(/i/r)^'^ scaling, we find frequent 
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out bursts which are triggered at small radii and which have slow rise times, in contrast 
to the observations. If the inner disk were missing due to a slow removal of matter onto 
the accretor, then the outburst would have to be triggered further out. This may not be a 
unique solution to the problem, and at least three obvious alternatives are apparent. (The 
problem of the quiescent X-ray flux would only be addressed by the third option.) 

(1) The constraint imposed by observed exponential decays a. — bQ{h/rY-^ is actually 
only apphcable to the matter which resides in the "hot state" part of the disk. Once 
the matter has made the transition to the low state, it does not participate in the flow 
dynamics. Therefore it is quite probable that a — 50(/i/r)^'^ is a scaling which applies only 
to ionized gas at T > 10^ K, and for neutral gas at T < 10^ K there may be another scahng 
which drops a abruptly to a value less than 50(/i/r)^-^. In this case, the material would 
evolve less in quiescence, and there would be less tendency for the matter to spread and 
trigger an inside-out outburst with its attendant slow rise. 

(2) We have not yet activated and tested the section of our code which provides for a 
variable outer disk radius (see Ichikawa & Osaki 1992). The quahtative effect this has on 
the evolution is to sweep up matter at large radii which shrinks the outer disk in quiescence. 
The contraction is caused by the accretion of low angular momentum material. This has 
the effect of enhancing the surface density near router and leaving the disk more prone to 
outbursts triggered at large radii (see Ichikawa & Osaki 1994). 

(3) There may exist a mechanism which prevents the cooling wave from propagating 
once it reaches some radius Tcrit- Since Mniin(r) = M(Einin) which provides a measure of the 
quiescent mass flow rate at a radius r scales as ~ r^'^, if rent were at ~ 10^ — 10^° cm, for 
example, this would be sufficient to leave a small, hot inner disk in quiescence that could 
give rise to the quiescent X-rays at a level ~ 10^° — 10^^ g s~^ seen in A0620-00 using a 
standard accretion efficiency e ~ 0.1. 
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Having made these concessions, we must return to our original point that the EUV/soft 
X-ray observations seem to demand some evaporation of the inner disk (unless the third 
option just discussed turns out to be true). Although we are uncertain of the physics 
which leads to evaporation, we should include at least an ad hoc prescription for removal of 
matter from the inner disk in any complete time dependent model. We therefore evacuate 
matter from the inner disk through the prescription dT,{r)/dt = — Eo/(l + [^/''"o]^)^) 
where Eq = 10'^{Mi/Mq)-^ g cm'^ 8"^ and ro = lOOrg = 2.95 x 10'^{Mi/Mq) cm (R. 
Narayan 1996, private communication). This gives an integrated mass removal Mevap = 
-J{d^/dt)27rrdr = 7r^or^J{l + [nn„er/ro]') = 2.7 x 10^^ g s'^ (1 + [ru,^,Jro]^y\M,/MQ). 
We must have some matter at each grid point for our code to work properly, therefore when 
E drops below 1 g cm~^ at a given radius we set it equal to 1 g cm~^. 

Figure 8 shows a series of light curves with and without evaporation. We use the 
expression given in the previous paragraph for Mevap, with the normalization multiplied by 
10"^. For the four panels shown, x — oo, A, 3, and 2. Hence in the top panel we present our 
standard model without evaporation. Outbursts arc triggered frequently — about every 5 
yr — at small radii in the disk. Most of the outbursts are fairly low amplitude, but since 
the 1975 outburst of A0620-00 made it a sustained ~ 50 Crab source in 3 — 6 keV X-rays as 
seen by Ariel 5 (Elvis et al. 1975), the fainter outbursts produced in our "standard model" 
would have been seen over the past ~ 20 yrs, had they occurred. About every 50 yr we get 
a major outburst in the models. The introduction of evaporation removes the inner parts of 
the disk, eliminates the minor outbursts, and increases the recurrence times for the major 
outbursts. Figure 9 shows the mass of the accretion disk accompanying the runs shown in 
Fig. 8. During a major outburst, a fraction ~ 0.2 of the disk is accreted, whereas during a 
minor outburst < 0.01 of the disk is accreted. 

Figure 10 shows expanded views of one major outburst taken from each of the four 
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panels in Fig. 8. The quiescent values of the accretion rate in third panel of Fig. 10 show 
the effective rate of the evaporation or mass loss from the disk. In the models of Narayan et 
al. (1996) and Esin, McClintock, & Narayan (1997), for instance, the efficiency of accretion 
is low in the quiescent state e ~ 10~^, while in outburst one still has e ~ 0.1. Therefore one 
would see a greater dynamic range in soft X-ray flux between quiescence and outburst than 
would be inferred directly from the rates of accretion given in Fig. 10. Also, our condition 
that E(r) > 1 g cm~^ reduces the effective rate of evaporation increasingly for smaller x 
values (i.e., higher evaporation rate): the effective evaporation rates in going from a; = 4 to 
X — 2 increase by a factor ~ 5 for each of the two factor of 10 increases in 10~^. 

In Fig. 10 we see that, as expected, the fact that triggering of outbursts cannot occur 
at small radii due to the truncation of the inner disk forces the rise times to shorten. Even 
for the highest values of Mevap; however, the rise times are ~ 1 month — longer than the 
value of several days commonly seen. For low values of the evaporation, only the quiescent 
surface density proflle is affected. The reason the outburst E(r) proflle is not strongly 
changed is that during outburst the surface density near the inner edge is so large that 
the local mass flow greatly exceeds the local evaporation rate. In looking at the decays 
of the outbursts in the last panel of Fig. 10, however, we see that for a high evaporation 
rate, the outburst itself is also affected: the decay is cut off so that the exponential decay is 
truncated. This occurs because, at some point during the inward movement of the cooling 
front, there is no longer any disk left into which the cooling wave may propagate. For 
instance, lor x — 2 we only observe about 2 decades of exponential decay in the value of M 
at the inner edge, which presumably powers the soft X-ray flux we observe. In systems like 
A0620-00, however, we observe about 3 decades of exponential decay (Tanaka & Shibazaki 
1996). The model with x — 2 corresponds to Mevap ~ 2 X 10^^ g s-^ with our standard 
assumed secondary star mass transfer rate value Mr ~ 2 x 10-*^^ g s"-*^. The rise time is 
~ 1 month. The rise time could be shortened presumably by increasing Mgvap further (i.e.. 



-19- 



decreasing x) , but that would reduce the dynamic range of the soft X-ray flux accompanying 
the exponential decay to an even more unacceptably small value. 

Figure 11 shows the evolution of E and Tmidpiane covering an outburst in the third panel 
of Figure 8. There are 10 time steps shown, with the first time step corresponding to just 
before the triggering of the outburst. The time increment 6t = 0.1 yr. The starting time 
step shows the truncated nature of the disk due to evaporation. The outburst triggers at 
logr(cm)~ 9.6. In the early stages of the outburst evaporation does not have any effect 
on the outburst, as the S(r) values are quite large for r small. In looking at time steps 
8, 9, and 10, however, we see that in the late stages of the outburst when E(r) becomes 
smaller, the evaporation swallows up the last of the inner, hot disk before the cooling front 
can propagate all the way to the inner edge. In time step 8 for instance the hot disk exists 
only in a ring from about 8.5 to 9.5 in terms of logr(cm). For time steps 9 and 10, the 
entire disk is back into the low state. This leads to the truncation of the full exponential 
decays evident in Figure 10. Figure 12 reveals the evolution in quiescence by showing four 
time steps which follow the last profile in Figure 11, here taking St = 10 yr. One can sec 
evaporation eating into the quiescent E(r) distribution, even as matter in the inner disk 
evolves viscously, and matter at the outer edge continues to arrive from the secondary star. 

In summary, evaporation appears not quite capable, at least for a system with 
parameters relevant for a system like A0620-00, of producing outbursts which occur at 
large enough radii to give rise times that are as fast as observed, and which simultaneously 
have ~ 2 — 3 decades of exponential decay. Since the models are so close to working, one 
could argue that, given the crude nature of the models, we could consider this a qualified 
success. On the other hand, evaporation unquestionably does help in the standard model 
as regards the recurrence times for outbursts. The effect or necessity of evaporation in the 
models may potentially be different for other model parameters, for instance with systems 
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for larger orbital periods and higher Mt values. HLMN find, for example, that evaporation 
is needed to obtain light curves in various energy bands that can account for observations 
of the start of an outburst seen in GRO J1655-40, a system with a much longer orbital 
period Porbitai = 2.6 d (Bailyn et al. 1995) and higher mass transfer rate Mt — 2 x 10^^ 
g than we have assumed in this work for our standard model (Porbitai — 7.75 h giving 
^outer = 1-5 X 10^^ cm and mass transfer rate My = 1.89 x 10^^ g s^^) taken to be 
representative of A0620-00. HLMN do not show light curves covering complete cycles of 
quiescence and outburst, so it is difficult to assess the full generahty of their conclusions. 
In particular, it may prove difficult for any disk instability model to reproduce the prompt 
decay in V which was seen in Hubble Space Telescope observations that occurred even as 
the 2 — 10 keV X-ray flux was increasing for the 1996 outburst in GRO J1655-40 (Hynes 
et al. 1996). The fact that the two microquasars do not show the classical fast-rise, 
exponential decay outburst may be associated with the fact that the BHs in these systems 
seem to be maximally rotating (Zhang, Cui, & Chen 1997). Our results do not support the 
strength of HLMN's contentions about the necessity of evaporation. Prehminary work using 
our canonical SS Cygni model (C93b; Cannizzo 1996c, 1998c) also flnds evaporation to be 
ineffectual in bringing about fast rise outbursts for models which are intrinsically prone to 
slow rise outbursts. Cannizzo (1998c) finds that, even in the limit where Mevap — ^ Mt, the 
outburst rise times as taken from light curves computed using Minner are still slower than 
the fast-rise outburst in SS Cygni observed by Mauche (1996) in the EUV. One could not 
of course have Mevap — Mt or the entire disk would evaporate. 



2.5. Strength of the Irradiation 

Tuchman et al. (1990) present a formalism for incorporating the effects of irradiation 
on the equilibrium structure of the disk determined by the vertical structure calculations. 
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We have included their formahsm into our time dependent model to calculate the local 
irradiation temperature and its effect in altering the S-curve. In this preliminary work we 
only activate the portion of the code which calculates the irradiation temperature. We leave 
for future work a computation which includes its feedback into the structure of the disk. 

Shakura & Sunyaev (1973) give the local irradiation flux as 

where C = Ad{h/r)/dlnr, n = 1 for neutron star accretors and n = 2 for black hole 
accretors. The parameter A is the X-ray albedo. The values of n advocated by Shakura & 
Sunyaev are based on the supposed preferred directions for the emission: isotropically from 
the neutron star surface (n — 1), and primarily vertically from the disk surface for black 
hole primaries (n = 2). One cannot claim to understand the emission geometry at small 
radii in the BHXBs — for instance there may exist a scattering medium above the disk 
which influences the value of n — therefore we shows results for both n — 1 and 2 as being 
representative of possible values. 

Figures 13 and 14 show the values of ratios of the locally deflned irradiation temperature 
Tirr = (Firr/cr)^/'^ to the local (viscous) effective temperatures for 600 d of evolution with 
St = 6 d spanning the rise and decay of an outburst in the standard model. We show 
n — 1 and n — 2, and assume A — 1. As in section 2.1, we adopt a logarithmic grid 
spacing so as to resolve better the inner portions of the disk. We take = 81 grid points. 
For the case n = 2 which is thought to be more representative of the BHXBs, the local 
irradiation temperature is never more than about 0.3 — 0.4 of the local effective temperature 
determined from equating viscous heating and radiative cooling. This occurs near the time 
of maximum light in the outburst. Van Paradijs (1996) presents revised criteria for the 
allowability of the limit cycle mechanism in the soft X-ray transients including neutron star 
and black hole systems. He concludes that irradiation is more effective for neutron star 
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accretors than for BHXBs at preventing the hmit cycle from operating by virtue of keeping 
the entire disk in the ionized state. In his Fig. 2 he indicates A0620-00 as lying a factor 
> 10 below the dividing line between steady and unsteady behavior. Therefore our finding 
of relatively weak irradiation for a model meant to represent A0620-00 is consistent with 
the fact that A0620-00 has exhibited two outbursts (in 1917 and 1975) and is seemingly 
not prone to irradiation induced effects. It is also consistent with HW and MW who found, 
by comparing theory with observation for A0620-00, that irradiation does not contribute 
overwhelmingly to the optical flux. 



3. DISCUSSION 



We have investigated a variety of phenomena associated with the limit cycle model as 
it may apply to the BHXBs. We find the following: 

(i) The grid spacing in which Ar oc r does not appear to be as convenient as the 
spacing Ar oc ^/r because the former overresolves the innermost part of the disk, resulting 
in long CPU times for execution. Low resolution runs with A?" ~ 20 — 40 grid points 
produce light curves which show artificial "steps" in the decay caused by the turn-off of 
individual grid points. In addition, the locally defined e— folding decays times are too slow 
by ~ 2 — 3, and the decay character is artificially slower-than-exponential. As N is made to 
increase, the steps becomes smaller and disappear, and the decay time scales asymptote to a 
common value, as expected in a numerical scheme which converges properly. We have only 
investigated properties associated with the decay. There may be other spurious systematic 
effects associated, for instance, with the long term light curves (i.e., the triggering radii 
and outburst rise times), or with the production of refiares or spikes in the light curve (e.g. 
K96) which would be smoothed out in a much higher resolution model. 
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(ii) The light curves for outbursts have slow rise times when triggered at radii < 10^ 
cm in the disk, and fast rise times when triggered at radii > 10^° cm. There is an abrupt 
change in the character of the rise at some radius intermediate between these two values. 
In this context, it is somewhat of a misnomer to refer to the two cases as "inside-out" and 
"outside-in" since, in both instances, the triggering radius is much less than the outer disk 
radius. The outburst light curves of interest in this work are typified by A0620-00 and show 
not only an exponential decay with ~ 30 — 40 day e— folding time scale, but also a fast rise 
(i.e., several days). These experiments in which we artificially trigger the outbursts at a 
specific radius show that, to produce outbursts with rise times as fast as those observed, 
the triggering radius must be at least as large as some value between about 10^ and 10^° 
cm, for A0620-00 parameters. 

(iii) The model of VW for the decay of the outburst, based on considerations involving 
a deviation from steady-state fiow in the inner part of the disk, is substantially borne out 
in detailed testing. Their basic premise is that the local fiow speed "accelerates" from 
the local value of aCs{h/r) to the local value of aCg within the hot part of the disk. VW 
conceive of the point at which Vr ~ aCs is attained as being at the inner edge of the cooling 
transition region, whereas we find it to be at the outer edge. This does not impact on VW's 
derivation. More importantly, we find that the value of the radius at which Vr changes sign, 
expressed in units of r(ihot)) asymptotes to ~ 0.4, whereas the fractional radius at which 
the local fiow rate equals the standard viscous fiow rate asymptotes to ~ 0.55. 

The work of CCL had as its underlying premise the assumption that the decay of the 
soft X-ray fiux in the outbursts seen in the BHXBs is due to a coohng wave. With that as a 
working hypothesis, CCL showed that cthot ~ 50(/i/r)^-^ is needed, and the e— folding decay 
time for dM/dt is ~ {)A{GMi/ aQ)c^^ where ~ 16 km s^^. An alternative starting point, 
advocated recently by King & Ritter (1998), is that the decays in the BHXB outbursts 
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are due to the gradual loss of disk material into the BH. The coohng wave is unable to 
start because of strong irradiation, and so the hot, viscous disk can only lose material from 
its inner edge. Cannizzo, Lee, & Goodman (1990) showed that, for a constant, dM/dt 
decreases in a purely viscous disk as ~ t'^-^. To get an exponential decay as observed, 
the kinematic viscosity coefficient v oc aTmidpiane must be constant, implying either an 
isothermal disk with a constant, or else a oc 1/Tmidpiane- Neither option seems physically 
reasonable. It is interesting that for both the cooling wave model and the viscous decay 
model, a constant cthot value does not work (Cannizzo 1998b). Mineshige, Yamasaki, & 
Ishizaka (1993) pointed out that purely viscous decays lead to power law light curves and 
that a cooling wave which extracts mass and angular momentum from the hot, viscous disk 
can produce the required exponential decay. They present a toy model in which dT^/dt 
oc —Tl/to is put in by hand, rather than determining a form for cthot that will lead to the 
exponential process self-consistently using a full time dependent model as in CCL. King & 
Ritter (1998) note the discrepancy between the X-ray to optical ratio in the computations 
of HW, and the observed value for A0620-00 in the 1975 outburst. The main problem is 
that HW's calculated peak M(rinner) value is ~ 3 x 10^^ g s~^ — a factor ~ 30 higher than 
observed. HW employed an early, experimental version of Cannizzo's time dependent disk 
code (from 1984), which may have led to numerical problems. For instance, they were 
forced to adopt an inner radius ~ 10 times too large. In this paper we use a much more 
modern version of the code and find M ~ 10^^ g s^^ as observed. 

(iv) In running long term models with the parameters given in CCL so as to cover 
complete cycles of quiescence and outburst, we find that outbursts tend to be triggered 
frequently at small radii, confirming the analytical estimates of Lasota et al. (1996). 
Small outbursts occur every ~ 5 yr, whereas major outbursts that could be observed with 
typical all-sky monitors on X-ray satellites occur every ~ 50 yr. The rise times for the 
outbursts are ~ 10 times longer than those commonly observed. When we include some 
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evaporation of matter from the inner edge so as to produce a truncated E(r) profile in 
quiescence, we find that the smaller outbursts are eliminated and the rise times associated 
with the major outbursts are shortened, but not enough to bring them into accord with 
observations. Therefore we cannot confirm the conclusions of HLMN. The parameters for 
their system of interest, the microquasar GRO J1655-40, are sufficiently different from 
A0620-00 that a direct comparison of our respective findings may not be justified. In our 
modeling we find that, if a given model for the quiescence is not able on its own to prevent 
inside-out outbursts from occurring, then the addition of evaporation as an outside agent 
will not influence the S(r) profile to the extent required. We note that the form for a 
arrived at by CCL is really only a constraint on a in the hot state (i.e., on the upper stable 
branch of the "S" curve). Were one to artificially lower a in quiescence by, for instance, 
reducing the proportionality constant ckq, it may be possible to force the matter to keep 
from evolving significantly at small radii and thereby prevent inside-out outbursts from 
triggering. Another possibility might be for the cooling front to "stall" at some radius 
~ 10^ — 10^° cm so as to maintain a permanent hot, inner disk. This would deplete the 
E(r) profile for small radii and produce quiescent X-rays at a level commensurate with that 
seen in A0620-00 using a standard efficiency e ~ 0.1. 

For our standard model plus evaporation at a level Afcvap ~ 10^"^ g s~^, wc produce 
outbursts with about the right spacing and amplitude as observed (trecur ~ 50 — 100 yr, 
outburst energy 5£' ~ 3 x 10^^ ergs, accreted mass SM ~ 5 x 10^^ g, and peak absolute 
visual magnitude My ^ +1), although the rise times are too slow. As noted earlier, Lasota 
et al. (1996) point out that the standard model for the limit cycle instability in the BHXBs 
produces outbursts which recur too frequently because the inner disk radius is quite small, 
and they argue for evaporation to truncate the inner disk. If evaporation at a level ~ 10^^ 
g s~^ is correct, then the efficiency of accretion in the low state must be e ~ 10~^ (Narayan 
et al. 1996) to account for the observation of McClintock et al. (1995) of A0620-00 in 



-26- 



quiescence. Following this one step further, because the efficiency in quiescence of the 
neutron star SXTs is still e ~ 0.1, there should exist a difference of ~ 10^ between the 
quiescent fluxes in BHXBs and neutron star binaries. This does not appear to be the case 
(Chen, Shrader, & Livio 1997, see their Fig. 14). Hence, while evaporation appears to be 
attractive in some ways, it may prove to be problematic in others. 

(v) The local values of the irradiation temperature are significantly less than the local 
values of the disk temperature in our standard model. King & Ritter (1998) argue quite 
the opposite, that irradiation is in fact strong during the outbursts of the BHXBs. They 
propose that the slow, exponential time scales for decay seen in A0620-00 and other similar 
systems are caused by the inabihty of the coohng front to form and shut off the fiow of 
matter onto the BH. We may estimate the strength of the irradiation during outburst in 
A0620-00 in several ways. Warner (1987) found an empirical relation between absolute 
visual magnitude My and orbital period Porbitai for dwarf novae at maximum light in an 
outburst. A similar relation was found by van Paradijs & McClintock (1994, hereafter 
vPMc) for outbursts in BHXBs and neutron star binaries. Their relation was between My 
and {Lx/ L-EddV^'^PoThitaX- The underlying cause for the Warner relation is thought to be 
the fact that the mass which can be stored up in quiescence in the limit cycle model is 
bounded from above by a critical value / 27rrdrllmax, and this value scales with the size of 
the system (C93b; Warner 1995; Cannizzo 1998a). There seems to be general agreement 
that the triggering mechanism for outbursts in both classes of systems is the limit cycle 
instability, therefore the "maximum mass" concept should apply to both. By comparing the 
My values between the Warner and vPMc relations, one sees that the SXTs are intrinsically 
brighter than the dwarf novae. VPMc argue that the form of their relation implies a strong 
X-ray heating effect which leads to the bulk of the optical fiux, through reprocessing, in 
the SXTs. In going between the dwarf novae and the neutron star SXTs, both of which 
presumably have ~ IM© accretors, the comparison at a given orbital period should be 
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direct and allow a pure estimate of the strength of the irradiation, because the sizes of 
the two systems and hence the surface areas of the accretion disks should be similar. For 
the BHXBs, however, one also needs to take into account the increased size of the disk, at 
given orbital period, due to the presence of a more massive central object. If we compare 
the Warner and vPMc relations at the approximate orbital period of A0620-00, we see that 

1 /3 

AMy ~ 3 mag. For Porbitai constant, the orbital separation scales as (if Mi » M2), 

so the radiating area of the disk which gives rise to the optical scales as M^^^ . In going 
from a ~ 0.7 — IM© WD to a ~ IMq BH, we increase the disk area and hence optical flux 
by ~ 7^/^ — 10^/^ ~ 4 — 5, which reduces the difference AMy between the Warner and 
vPMc relations for an A0620-00 type system by ~ 2 mag. Therefore AMy (corrected) ~ 1 
mag. This comparison is admittedly crude because there may also exist a correction due to 
an effective temperature difference (J. van Paradijs 1997, private communication). We can 
obtain a more direct estimate simply by comparing our computed My for A0620-00 with 
that observed in the 1975 outburst. The shortfall in optical flux in the models which ignore 
irradiation must be due to the reprocessing of X-ray flux. Therefore by determining the 
deficit, we learn how much of the optical flux comes from reprocessing, and how much is 
due to viscous dissipation. This exercise has already been carried out by previous accretion 
disk hmit cycle modelers specifically for A0620-00 using time dependent models which 
did not include the effects of irradiation. HW modeled the B band flux and determined 
AMb ~ 1 mag at maximum light (see their Fig. 3), and MW modeled the V band flux 
and found no difference between their computed value and the observed value at maximum 
light (see their Fig. 10). In our model we find, for a face-on disk. My ~ -|-1 at maximum 
light, whereas the observed value, corrected for extinction, is -1-0.7 (vPMc). If we increase 
the inclination of our system to ~ 70° (Haswell et al. 1993), our system becomes fainter 
and we add ~ 1 mag to our My (peak) value, thereby giving AMy ~ 1 between theory and 
observation, in line with the other estimates. Finally, in comparing our nonirradiated and 
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irradiated models, we find at the peak of the outburst that Tirradiation/T'effective ~ 0.3 — 0.4 
for n = 2 thought to be applicable to the BHXBs. VPMc show that, for the temperature 
regime of relevance for the SXTs, the optical flux scales roughly as Tjj. Therefore were we 
to self-consistently include the effects of irradiation in computing My, we would produce 
more optical flux by a factor ~ 1.3^ — 1.4? ~ 2, which would subtract ~ 1 mag from the 
inclination corrected value My ^ 2, bringing it back into accord with the observed value. 
In summary, it seems that the contribution of the reprocessed X-ray flux to the optical flux 
at peak hght in A0620-00 is ~ 1 mag, meaning its relative contribution is ~ 50 — 70%. This 
can be produced in the time dependent disk models by a slight increase of the effective 
temperature in the inner part of the hot, viscous accretion disk. Were the entire disk to 
be strongly irradiated out to ~ 0.7 of the Roche lobe of the primary so as to prevent the 
coohng wave from forming, as advocated by King & Ritter (1998), the optical flux would 
be overproduced with respect to what is observed in A0620-00. 

We may quantify this last statement using the subroutine of our time dependent code 
which calculates the V band flux. If the entire disk is irradiated so that the cooling front 
is unable to form and propagate, that means the effective temperature must be at least 
~ 9000 K. The irradiation flux falls off as r~^, so the radial distribution of the effective 
temperature, if it is dominated by irradiation flux, is r~^/^. Inputting the law Tcflf(r) = 9000 
K (r/ 

^outer) ) where router — 1.5 X 10 cm, we flnd My — +0.45. For a moderate 
inclination, this becomes ~ 1.5. If the entire exponential decay observed in V during the 
1975 outburst of A0620-00 were irradiation dominated, this is the faintest the disk could 
become. At the end of the exponential decay, however, one saw my ~ 14 — 15 (Lloyd 
et al. 1977), corresponding to an absolute V magnitude My ~ 3 — 4 (vPMc) — fainter 
than seemingly allowed by the irradiation model. As a test of the zero point of the V 
calibration, we input T^s = Tq = 5770 K (constant with radius) and router = R& = 7 x lO-*^*^ 
cm, and flnd My = +4.86. The accepted solar value is My = +4.83 (Allen 1976). King 
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& Ritter (1998) point out that the constraints on the functionahty of a. derived from 
arguments hke those presented in Cannizzo et al. (1990) do not apply for disks in which 
the midplane temperature is determined not by viscous heating but rather by irradiation. 
For the midplane temperature to be strongly influenced, a rather large irradiation is needed 
(Tuchman et al. 1990). Our estimate that ~ 10^ K at the outer disk edge may be 
unrealistically low. During the time in which the optical flux decreased by ~ 2 — 3 mag in 
the 1975 outburst of A0620-00, the X-ray flux fell by a factor ~ 10^. It is not clear how 
the strongly irradiated midplane would be able to maintain a constant temperature Tmid 
during this period of declining irradiation so as to lead to the kinematic viscosity coefficient 
V being constant with time in the outer disk. The constancy in v is needed in order to 
produce an exponential decay of the visual and X-ray fluxes in the model of King & Ritter 
(Mineshige et al. 1993). 

In this work we have chosen to concentrate on the most basic features of the BHXB 
outbursts: the recurrence time scale, the ~ 30 d exponential decay, and the < 3 d fast 
rise. We do not directly address the issue of the secondary maxima or refiares which are 
frequently seen. Based on some of our findings, however, we may discuss previous work 
on the origin of the refiares. These are resurgences in the X-ray flux which are also seen 
to some extent in the optical, during which time the flux increases by a factor ~ 2 before 
resuming its decline. The decay properties after the reflare are much as before, so it is as 
if the light curve were displaced upward by a factor ~ 2 beginning at a point ~ 60 — 70 
d after the outburst starts. Several works have appeared concerning the refiares (Chen, 
Livio, & Gehrels 1993, hereafter CLG; Augusteijn, Kuulkers, & Shaham 1993; Mineshige 
1994; K96; Kato, Mineshige, & Hirata 1995, hereafter KMH; IMK; Kuulkers, Howell, & van 
Paradijs 1996, hereafter KHvP; for a review see Wheeler 1998). One feature which some 
of the theories of the reflare share is that they involve irradiation, either of the secondary 
star so as to induce additional mass overflow which powers an increase in the luminosity by 
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virtue of the extra material added to the disk (CLG; Augusteijn et al. 1993), or irradiation 
of the disk so as to reflect the cooling wave at one or two grid points as a heating front and 
momentarily rebrighten part of the disk (K96). 

In the context of the theories which rely on irradiation, it is interesting that the reflares 
seem to be uniquely associated with the BHXBs (Wheeler 1998). Reflares are not seen 
in the neutron star SXTs such as Cen X-4 and Aql X-1. According to our estimates of 
the strength of irradiation, and that of previous workers (HW; MW; van Paradijs 1996), 
however, irradiation is weak in the BHXBs. If reflares arc unique to BHXBs, and if 
irradiation is weaker in BHXBs than in neutron star SXTs, it seems rather unlikely that 
irradiation can be the BHXB specific parameter which causes reflares to occur only in the 
BHXBs. This makes it difficult to accept any of the irradiation-based models. Another 
consideration is that, whatever the physical mechanism is, it must be such that it leads to 
a delay time between onset of the main outburst and secondary reflare that scales with 
orbital period. The fact that the reflare in GRO J0422-I-32 occurred ~ 45 d after primary 
maximum, whereas in A0620-00 it occurred ~ 55 d after primary maximum, led CLG to 
predict an orbital period less than 7 hr for this system. The period was later found to be 
5.1 hr (Ghevaher & Ilovaisky 1994, Orosz & Bailyn 1995). 

What unique property of BHXBs — something which would not be true for the neutron 
star SXTs — exists that could account for the reflares? This has already been answered by 
many authors — it is the extreme value of the mass ratio q = M2/M1 (Bailyn 1992; KMH; 
IMK; KHvP). This ratio is much smaller for BHXBs than for neutron star SXTs, and this 
can lead to resonances occurring within ~ 0.7 of the Roche lobe of the primary wherein 
the disk is constrained by tidal truncation to exist. One observes superhumps (SHs) during 
superoutbursts (SOs) in the SU UMa class of dwarf novae (Warner 1995). The SHs are 
thought to be caused by the fact that the 3:1 commensurability resonance with the binary 
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orbital period lies within the accretion disk for q < 0.25 (Whitehurst 1988, 1994, 1995), but 
lies beyond the point of tidal truncation for q > 0.25. In fact, SHs are only seen in short 
orbital period dwarf novae for which q < 0.25. Whitehurst showed that when the condition 
q < 0.25 is fulfilled, the outer disk acquires an eUiptical shape which precesses in the 
corotating frame of the binary, and this precession causes the distance through which the 
mass stream must fall before striking the outer disk to vary, thereby modulating the light 
curve on a time scale equal to the orbital period plus a few percent. Osaki (1989) combined 
the disk instability model with a proposed "tidal instability" which he envisioned to occur 
whenever Whitehurst's condition for superhumps becomes satisfied. Osaki's thermal-tidal 
(TT) instability allows for superoutbursts in the SU UMa stars by imagining that, whenever 
the outer edge of the disk expands beyond the 3:1 resonance in the disk, the increased 
tidal torque on the disk enforces a contraction of router by ~ 10 — 20 %. The mass that 
was present in the disk at the start of the outburst is confined to a smaller volume, and 
the increase in Xl(router) relative to Sinin(?^outer) produces an "over-filled" disk that prevents 
the cooling wave from starting until a sizable fraction of the disk matter has accreted onto 
the central object. The motivation for the contraction of the outer disk was not obvious, 
based on Whitehurst's time dependent calculations, and led to criticisms (Whitehurst 1994; 
Cannizzo 1996a). Recent work by Murray (1997) which utihzes a more advanced SPH 
model than was available to Whitehurst, however, does in fact show a contraction of the 
outer radius which is coincident with the growth into the nonlinear regime of the elliptical 
disk induced by the 3:1 resonance. This gives credence to Osaki's model. 

What does this have to do with the refiares in the BHXB outburst light curves? A 
potentially key observation is that, in the BHXB GRO J0422-I-32, SHs appeared in the 
outburst at about the same time as the reflare, thereby arguing for a direct connection 
between the nonlinear outcome of the 3:1 instability — namely the formation of an eccentric 
disk and concomitant SHs — and the refiare (KMH, KHvP). The crucial difference which 
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the global effect the disk contraction has in going from accretion disks in dwarf novae to 
those in the BHXBs has to do with the growth time to reach the nonlinear regime of the 
tidal instability. Several workers have suggested a direct connection between the SHs and 
the reflares (KMH, IMK, KHvP), namely, that it is the contraction of the disk in Osaki's 
TT model that leads to the reflares. IMK argue that the 3:1 resonance growth time varies as 
^grow ~ l/(^['"outer]), wherc fl^Touter) IS the Kepleriau frequency at the outer disk edge. They 
point out that igrow is about 30 times larger for the BHXBs than for the SU UMa stars. 
For normal outbursts which turn into SOs in the SU UMa stars by virtue of the expansion 
of router beyond the 3:1 resonance radius, the time tgrow(SU UMa)~ 1 — 3 d, so that one 
sees in the light curves first a rise to maximum associated with the normal outburst, then a 
brief dip, and then a rebrightening to maximum as the forced contraction of the accretion 
disk becomes active and the disk gas gets confined to a smaller radius. This constitutes the 
formal start of the SO, marked subsequently by a long term slower-than-exponential decline 
that typifies decays from maximum light in which the cooling front is unable to propagate. 
If tgrow ~ 100 d for the BHXBs, this would have important consequences for the interaction 
of the contraction phase with the inward propagation of the cooling front. This time scale 
is about three e— foldings for the time to deplete the mass of the hot part of the accretion 
disk, or one e— folding time \dlnrf/dt\~^ for the coohng front to move inward (CCL), so 
the cooling front has already had ample time to travel inward a considerable distance from 
^outer as defined by the tidal truncation, and subsequently redefined by the increased degree 
of tidal truncation due to the 3:1 resonance- induced elliptical disk. There is no opportunity 
for the cooling front to be prevented from forming as in the SU UMa SOs, and yet the gas 
in the accretion disk must be confined in some manner to a smaller value. It is difficult to 
estimate the degree of this contraction without doing a detailed calculation, considering 
especially the fact that the contraction occurs primarily in the outer, cold portion of the 
disk. One suspects that the brief period of augmentation in the local E(r) values would 
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produce a reflection of the cooling wave as a heating wave, which would propagate out some 
distance and then reflect inward again as a cooling front. This is found in time dependent 
calculations (IMK). The amplitude of the reflares observed in the BHXB outburst light 
curves are about a factor of 2 with respect to the X-ray flux observed just prior to the start 
of the reflare. This would not require a large reflection in order to be produced. The reflares 
caused by the reflection of cooling waves into heating waves seen in some of the light curves 
in CCL showed increases of ~ 10 in flux. (Their cause is thought to be numerical, however, 
and not physical.) In summary, the fast growth time for the 3:1 instabihty in the SU UMa 
stars leads to an immediate contraction in router and consequent over-flUed disk which can 
only decay viscously (since the cooling front cannot propagate), whereas the slow growth 
time for the 3:1 instability in the BHXBs means that the contraction, when it does occur, is 
too late to produce an over-flUed disk. Instead, one has only a minor re-adjustment in E(r) 
which produces a reflected transition front and resultant reflare. 

IMK attempted to calculate models using the TT instabihty for BHXBs, but various 
deficiencies in their work led to outbursts which did not resemble those observed. They 
use a low resolution logarithmic spacing ~ 20 — 40 to cover a dynamic range in disk 
radii router/^nner > 10^, they use a parametrization of a oc r°'^ (C94) which was criticized 
as not being valid for BHXB parameters (CCL), and their model contains detailed ad hoc 
prescriptions for the functional form of the increase in the tidal torque associated with the 
TT instability (both a and the torque proportionality constant each have ~ 6 adjustable 
parameters) . The hght curve they present covering an outburst with the TT induced reflare 
shows a factor ~ 10^ increase in Lx in the reflare, rather than the observed factor ~ 2 
increase. The decay characteristics are also not as observed. Nevertheless their model 
represents a flrst step toward a complete theory of outburst plus reflare in a time dependent 
model. Work is currently underway to activate and test the portion of our code which 
allows for a variable outer disk radius so as to test the global manifestation of Osaki's TT 



instability in the BHXBs. 
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4. CONCLUSION 

We have investigated several aspects of the accretion disk hmit cycle mechanism as 
it applies to BHXBs. Our main impetus is to extend work reported in CCL by running 
more extensive time dependent models, cover complete cycles of outburst and quiescence, 
and look in detail at several aspects of the model not studied by CCL. We find that the 
explanation for the cooling front speed put forth by VW seems to account for the deviation 
from steady state flow which we observe. In running complete cycles, we flnd a strong 
tendency using the CCL a form for frequent outbursts to occur which begin at small 
radii in the disk and have slow rise times, in contrast to the observations. Evaporation 
of material from the inner disk does not shift the surface density distribution enough 
to decrease the rise times sufficiently, although it does increase the recurrence times for 
outbursts to an acceptable value (Lasota et al. 1996). A separate issue is the necessity of 
having some mechanism distinct from the standard limit cycle operating in order to supply 
material to accrete onto the central object in quiescence so as to provide substantial X-ray 
and EUV radiation. Finally, for our standard model which is meant to represent a system 
like A0620-00, irradiation of the disk does not appear to be strong enough to change the 
structure significantly. This may pose a challenge for workers that seek to invoke irradiation 
to prevent the cooling front from propagating. Also, the apparent weakness of irradiation 
in the BHXBs relative to the neutron star binaries may present problems for theories of the 
refiare which are based on irradiation. 

This paper is not intended to be the final word in this subject, but should only serve 
to pose additional questions which must be looked into. There is clearly much more work 
to be done. In spite of our efforts we have failed to produce a successful "standard model" 
for A0620-00 which reproduces not only the ~ 30 d e— folding time exponential decay and 
~ 60 yr recurrence time scale for outbursts but also the fast rise — ~ 3 — 5 d. We have 
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purposely avoided introducing modifications to the a form utilized by CCL in order to 
provide continuity with that work. It now appears, however, that it may prove necessary to 
reduce a in the low state below that which would naturally occur using the CCL scaling. 
This could be carried out for instance by reducing ao. There may be other ways to get the 
fast rise, such as having a stalled cooling front which would deplete E(r) for small r. 

In this work we have only activated the portion of the our code which calculates the 
irradiation temperature based on the central flux from the disk and the local shape, as 
explained in Tuchman et al. (1990). The next step is to activate and test the portion 
which also "corrects" the local disk temperature (both effective and midplane) based on the 
irradiation temperature, which would in turn affect the global evolution of the disk (as in 
K96). Since the irradiation temperature does not even approach the effective temperature 
for parameters relevant to A0620-00, the kinematic effect on the disk would be insignificant. 
For the neutron star systems, however, irradiation will most likely play an important role. 

We have also chosen not to activate the part of the code which treats the variation of 
outer disk radius (Smak 1984). This may affect the nature of the hght curves accompanying 
the outbursts, in particular the rise times, by virtue of the mass swept up at large radii by 
the arrival of low angular momentum material from the secondary star. 

We acknowledge useful discussions with W. Chen, J.-M. Hameury, A. King, E. 
Kuulkcrs, J. -P. Lasota, C. Mauchc, J. McClintock R. Narayan, J. van Paradijs, E. Vishniac, 
B. Warner, and C. Wheeler. JKC was supported through the long-term scientist program 
under the Universities Space Research Association (USRA contract NAS5-32484) in the 
Laboratory for High Energy Astrophysics at Goddard Space Flight Center. 
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FIGURE CAPTIONS 

Figure 1. The decay from maximum light for a disk initially entirely in the high state, 
but with Xl(r) outer < Smin SO as to ensure the cooling front starting. We show the light 
curves in absolute visual magnitude M„ {top panel) and the locally defined decay time 
constant in days/magnitude {bottom panel). The six curves are for N = 21, 41, 81, 161, 
321, and 501. The last curve was started from an initial condition with the disk in the low 
state, but with a shght excess in E at 10^ (cm) to trigger the instabihty. 

Figure 2. The variation in local flow rate and local disk thickness during a decay. The 
plot shows Vr/aCg versus r {top panel) and h/r versus r {bottom panel). We show eight 
time steps during the decay spaced 10 days apart as the cooling front moves inward. The 
x's indicate individual grid points. The dashed line in the upper panel shows the locus of 
grid points ihot which lie at the outermost point of the inner, hot disk. The dashed line in 
the second panel shows the values oi h/r associated with Smin, and the dotted hne gives 
conditions at Smax- 

Figure 3. The radii at which the local flow speed (i) goes to zero and (ii) equals the 
standard value expected in a steady disk aCs{h/r). Both are given in units of the radius of 
the hot portion of the disk r(ihot)- The eight time steps shown correspond to those depicted 
in Fig. 2. 

Figure 4. The local time constant associated with the decay of the optical flux for 
eleven values of n, where a — ao{h/r)". The value of n is 1.50, 1.52,..., and 1.70. To keep 
the value of My within a common range, the proportionality constant logao is varied in 
steps of 0.04 from 1.7 to 2.1. 

Figure 5. The variation of properties associated with the transition front during 
an outburst which is triggered at a small radius. We show the number of grid points 
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with the front A^/, the fractional width of the front w/r, and the ratios \Jh/r and h/r 
evaluated in the outermost radial zone within the hot state. Within the rise of the outburst 
w/r ~ (/i/r)°'^^, whereas during the decay of the w/r ~ {h/r)^^^. 

Figure 6. The dependence of the outburst profile on the radial location of triggering. 
The panels show the initial low state configuration {top panel), and the resulting light curves 
{bottom panel). The three small Gaussians in the top panel are representative of the narrow 
surface density enhancements which are added to the broad, underlying configuration, in 
order to trigger the outburst. The functional forms of the narrow and broad Gaussians are 
given in the text. 

Figure 7. The variation of the surface density with radius and time for inside-out 
outbursts and outside-in outbursts. We show a model run for log reenter (cm) = 9 and for 
log reenter (cm) = 10. The dashed curve gives the initial state, which is entirely along the lower 
branch of the S curve, and the dotted lines correspond to the local maxima and minima in 
E from the steady state scalings. We show 30 d of evolution, with a 5 d spacing between 
each curve. In the upper panel there are six sharp spikes evident between logr(cm)= 10 
and 11 showing the outward movement of the heating wave at t{d)— 5, 10, 15, 20, 25, and 
30. In the lower panel there are also six spikes in a similar radial range showing the outward 
movement, as well as five spikes between logr(cm)= 9 and 10 showing the movement of the 
inner edge of the heating front at t{d)— 5, 10, 15, 20, and 25. For t{d)— 30, the hot region 
extends all the way down to logr(cm)= 7. 

Figure 8. A series of outbursts spanning several cycles of quiescence and outburst. The 
t/— axis shows the log of the mass accretion rate from the innermost grid point in g s~^, and 
the X— axis gives the time in yrs. The top panel represents the standard model without 
evaporation. For the next three panels we remove material from the inner edge of the disk 
at the rate given in the text, with the normalization constant given in the text multiplied 
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by 10 ^, where x — 4, 3, and 2 {top to bottom). 

Figure 9. The time history of the mass of the accretion disk accompanying the four 
hght curves shown in Fig. 8. 

Figure 10. One major outburst taken from each of the four panels in Fig. 8. The 
values of x are given beside each curve. Smaller x values correspond to larger evaporation. 

Figure 11. The evolution of E and T accompanying an outburst in the third panel of 
Fig. 8. Ten time steps are shown, with the first corresponding to just before an outburst 
begins. There is a spacing of 0.1 yr between curves. The high temperature at the inner disk 
edge is an artifact of our imposition that E(r) > 1 g cm~^, which prevents the inner edge 
from cooling completely. 

Figure 12. The evolution of E and T continuing from Fig. 11, showing the quiescent 
period following the outburst. The first time step corresponds to the last step in Fig. 11. 
Four time steps are shown, and the spacing is 10 yr. 

Figure 13. The ratio of the local irradiation temperature to the local effective 
temperature for 600 days of evolution, shown at 6 d intervals, for n = 1. The drop-off at 
large radii moves inward as time progresses. 

Figure 14. The ratio of the local irradiation temperature to the local effective 
temperature for 600 days of evolution, shown at 6 d intervals, for n = 2. The drop-off at 
large radii moves inward as time progresses. 
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